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A remark on "On the Formation and Progenitor of PSR J0737-3039: New Constraints 
on the Supernova Explosion Forming Pulsar B" by Willems et al., [JjJ 

Tsvi Piran & Nir J. ShaviAQ 

Racah Institute of Physics, Hebrew University, Jerusalem 91904, Israel 

Shortly after the discovery of the binary system PSR J0737-3039 ^3, ll| we 13, 1^ predicted that it 
has a small proper motion and suggested that this implies that the progenitor mass of the younger 
pulsar (B) must be around 1.45 Mq. This is in contradiction to standard evolutionary scenarios 
that suggest a He-star progenitor with M > 2.1 — 2.3 Mq 0,0 and in sharp contrast to a velocity of 
~ 100km/ s predicted by these models 0- It requires a new previously unseen type of gravitational 
collapse. Our (generally ignored) prediction of a low peculiar motion was confirmed later by pulsar 
timing observations |g] that found v± < 30 km/s. Recently Willems et al. Q put forward an 
elaborate simulation suggesting that in spite of these observations, the most likely natal kick is large 
(50-170 km/s) and the corresponding most likely progenitor mass is large as well. We refute this 
claim and show that Willems et al. ||l| implicitly assume the result that they obtain. 



PACS numbers: 97.60.Gb,97.80.-d,97.60.Bw,97.60.- 



FolUowing the discovery of the binary system PSR 
J0737-3039 0, H we H Q reahzed that its location, 
^ 50 pc from the Galactic plane, implies that the velocity 
of the system is very small. This, in turn, poses strong 
limits on the natal kick velocity of the system and on the 
progenitor mass of the younger pulsar (B). The kinemat- 
ically favored progenitor mass is around 1.45 Mq. This 
low mass requires a new previously unseen gravitational 
collapse mode. Our prediction that the system has a 
small proper motion was in contradiction to the predic- 
tions of the standard model 0, 0| that assumed that the 
progenitor was a He-star with M > 2.1 — 2.3 Mq. This 
would have implied a typical proper motion of more than 
100 km/s. 

Our prediction of a very small proper motion was in- 
deed verified when pulsar timing set an upper limit of 
30 km/s on the proper motion of this system [l^l . In spite 
of that, Willems et al. p| presented elaborate simulations 
suggesting that the most likely natal kick is large (50-170 
km/s) as is the most likely progenitor mass. They recon- 
cile these suggestions with the limit on the proper motion 
by assuming that the system is moving with a very large 
velocity almost exactly towards us. Willems et al. 1'] in- 
sinuate that the difference between our conclusions and 
theirs arises because of the "higher accuracy" of their 
calculations that involve a three dimensional tracking of 
the galactic motion of the system whereas we @ consider 
only the vertical motion. We refute this claim and show 
that Willems et al. P| implicitly assume the result that 
they obtain. 

We begin with a brief simplified outline of our reason- 
ing. For brevity we skip important but non essential de- 
tails 0. If a significant mass is ejected during the second 
supernova collapse the binary system will be disrupted 
unless the second neutron star gets a large (comparable 
to its orbital velocity) natal kick in the right direction. 
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Getting the post SN system into an almost circular or- 
bit (as observed in J0737-3039) requires additionally a 
specially tuned direction. The kick given to the newborn 
pulsar will also be reflected by a large CM velocity, vcmo, 
of the binary system, which will be of order of a third of 
the orbital motion of the resulting system (this orbital 
velocity is - 600km/s for PSR J0737-3039). Determi- 
nation of vcMo, the CM velocity right after the second 
supernova can therefore set an upper limit on the ejected 
mass and from there on the progenitor's mass. 

It is impossible, of course, to measure vcmo- But if 
we can estimate the current CM velocity, vcm we would 
have a first approximation for vcmo- One can expect 
that in a typical galactic motion the magnitude of the 
planar components of peculiar motion will not change 
much 11]. To estimate vqai we used, in 2004, the lo- 
cation of the system Zobs ~ 50pc from the plane of the 
Galaxy. Stars move in a periodic motion in the vertical 
direction. For small vertical oscillations, the potential 
of the Galaxy is harmonic with a vertical orbital period, 
Pz « 50 Myr. The typical velocity for an object at Zobs 
is Vz ~ 2iTZobs/ Pz- Zobs ~ 50pc implies that the expecta- 
tion value of the vertical velocity is of the order of 6km/s. 
This reasoning has lead us to predict that both v_\_, the 
proper velocity of the system that was not known at that 
time, and vcm will be small. 

By now, this stage is not necessary. We know that 
f_L < 30km/s |3- This is a much stronger result than 
what we had in 2004. First the statistical uncertainty 
involved with the possibility that the system is just acci- 
dentally in the galactic plane disappears. Moreover, the 
vertical motion constrained just one component of the 
velocity whereas the present limit of v±^ constrains two. 

However, we still need to know the full magnitude of 
Vcm- Using Bayes' theorem we write: 

P{vcm\vi_ < 30km/s) cx P{vj_ < 3Qkm/s\vcM)P{vcM)- 

(1) 

Simple geometric arguments allow us to estimate the first 
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factor on the rhs: 

P{v± < 30km/s|t;cM) = 1 - cos P° '""'^" ^ . (2) 

V VCM ) 

If we are careful not to inflict our theoretical bias on the 
resulting distribution we should not have any prior on 
the probability distribution of vcm- In this case, as ex- 
pected a small observed v± implies a small vcm- We find 
VCM < 150, 100 km/s at 1% and 5% levels respectively. 
The most likely value is of course much smaller. This ba- 
sically leads to our 0, conclusion that vcmo was small 
from which it follows that the corresponding natal kick 
was small and finally that the progenitor's mass must 
have been small as well. In fact our best bet progenitor 
mass 1.45 Mq corresponds to ejection of 0.2 Mq (the ex- 
pected minimal mass loss in form of neutrinos) . Both the 
observed eccentricity and the limit on v± are compatible 
with such a collapse with no kick velocity at all. 

In order to estimate the effects of three dimensional 
Galactic motion Willems et al. Q integrate back- 
wards the orbital motion of the system from its present 
place towards its birth place. As such system are typi- 
cally born within the Galactic plane they identify possi- 
ble birth places as the intersection of the orbit with the 
galactic plane. Once the system's velocity in its birth 
place Vcmo is known, they proceed to estimate the natal 
kick and the progenitor's mass. Repeating many times 
they obtain a distribution of velocities and masses. 

This procedure, which seems reasonable at first glance, 
requires a knowledge of the current velocity and in partic- 
ular of the present radial component of the CM velocity. 
However, this is exactly what we are trying to infer from 



the data. Willems et al. Ill realize this and con- 
sider several distributions of the radial velocity. These 
distribution are either ad hoc (a uniform velocity distri- 
bution from -1500 km/s to -1-1500 km/s) or derived from 
a population synthesis model (Gaussian with a velocity 
dispersion 60 km/s - 200 km/s). 

In the language of our analysis these assumed distribu- 
tions amount to priors on P{vcm) that appears in Eq. 
n Clearly with a suitable prior one can get any result 
one wants for P{vcai\v± < 30 km/s) regardless of how 
small is the measured proper velocity and if vc m is large 
both the kick velocity and the progenitor's mass are large. 
Thus, Willems et al. simply assume (via this prior) 
the result that they obtain. 

There is not much to say about the assumption of an 
ad hoc uniform velocity distribution. It yields typically 
large velocities (less than a 10% chance of a velocity less 
than 100 km/s) forcing the distribution to favor large 
natal kicks and large progenitor masses. The distribu- 
tion derived from the population synthesis is more subtle 
as it presumably follows a valid theoretical justification. 
However, the population synthesis calculations are based 
on the current standard model for binary NS formation, 
and are hence implicitly biased towards large progeni- 
tor masses. Namely, all that the authors have done was 
to assume a large progenitor mass, simulate forwards in 
time to find that it implies a large radial velocity, that 
must be well aligned with the line of sight towards the 
observer. They then feed this velocity as input and inte- 
grate backwards in time and find, to ones great surprise, 
a large progenitors mass. At least the calculations are 
consistent. 
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